(N 

On 
ON 



Wave functions and their use in spectroscopy and phenomenology 

Thomas A. DeGrand a and Matthew W. Hecht b 

a Department of Physics, University of Colorado 
Boulder, Colorado 80309 USA 

b National Center for Atmospheric Research 
P. O. Box 3000, Boulder, Colorado 80307 USA 

We describe the calculation of Coulomb gauge wave functions for light quark systems, and their use as inter- 
polating fields for excited state spectroscopy. 
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1. WAVE FUNCTIONS 

We have been using wave functions recon- 
structed from lattice Monte Carlo simulations of 
QCD to calculate the masses of orbital excita- 
tions. We have also been trying to use wave func- 
tions as a tool for hadronic phenomenology. Our 
work has been presented in a series of papers El . 
We present a synopsis of it here. We will also 
comment on issues relevant to excited state spec- 
troscopy which came up during discussions during 
the conference. 



2. SPECTROSCOPY OF ORBITAL EX- 
CITATIONS 

Lattice QCD is unique among the subfields of 
physics in focussing only on ground state spec- 
troscopy. Some P-wave states' masses are regu- 
larly measured in staggered simulations because 
they are the odd parity partners of "ordinary" 
states: the a\ and p are examples of such pairs. In 
nonrelativistic QCD, Lepage and Thackerg] have 
computed the masses of \C a-ud Xb states (with- 
out including spin effects). The APE collabora- 
tion []3| measured masses of some P-wave mesons 
in quenched simulations at 6/g 2 = ft = 5.7, 
but had difficulty in continuing their program to 
higher (3 [Q. Recently, the Fermilab group has 
presented a calculation of the IP-IS splitting in 
charmonium, which they use to fix the strong cou- 
pling constant ||. 

We construct orbitally excited states by us- 
ing interpolating fields which couple only to a 
specific angular momentum eigenstate, which are 



projected onto zero momentum and which are of 
large spatial extent to maximize overlap with the 
state. 

At the t ^ end of the correlation function 
we use an operator which depends on the relative 
separation of the quarks, a "wave function" m . 
The wave function ipci^) of a hadron H in a gauge 
G is defined as 



i>G(r) = ^2(H\q(x)q(x 



r)|0) 



(1) 



where q(x) and q(y) are quantum mechanical op- 
erators which create a quark and an antiquark at 
locations x and y. (We have suppressed Dirac 
and color indices.) Our correlation function is 
constructed from convolutions of quark and anti- 
quark propagators G(x, y) 

C(f,t) = ^(ft.fcWjri.O;*,*) 

G q (y 2 ,0;x + r 7 t) (2) 

where ^(yi, j/b) is the t — operator. At large t 
if the mass of the hadron is run, then 



C(r,t) ~ exp(-m H t)ip G (r) 



(3) 



and so by plotting C(r, t) as a function of r we 
can reconstruct the wave function up to an overall 
constant. We measure the mass of a state by 
convoluting C(r, t) with some test function which 
further projects out the desired state: 

C(t) = X>Lrf(OC(?,t). (4) 

r 

At t = we take an operator which is separable 



in the coordinates of the quarks. For a meson we 
use 



^{x 1 ,x 2 ) = (j>i(xi)4> 2 (x2)- 



(5) 



In order to couple to orbital excitations we take 
4>i to be an S-wave and 4>2 to be some orbitally 
excited state with angular momentum I, centered 
around some specified coordinate. This state is a 
linear superposition of a p = L = I orbital exci- 
tation and a state whose center of mass momen- 
tum is nonzero. (This is the familiar "translation 
mode" of a shell model state.) Convoluting quark 
propagators as in Eqn. pi removes the p^O state 
and gives us the wave function of the p = L = I 
state. 

Our trial states 4>(x) were chosen to be Gaus- 
sians times an appropriate spherical harmonic. 
We used Coulomb gauge. We discussed with D. 
Richards whether P-wave spectroscopy could be 
done with other sources and concluded that one 
could use Wuppertal sources (see Ref. M): the 



S-wave 6 could be the inverse of D 



and 



the P-wave (f> could be cf>p = Di<&s- This would 
eliminate any need to gauge fix or perform con- 
volutions. We have to say, however, that we pre- 
fer the possibility of reconstructing the sink wave 
function at the end of the calculations; we have 
seen too many QCD simulations where one dis- 
covered at the end of a lot of data collection that 
one's operators were not as effective at producing 
the desired state as one would have wished. 

There were several things that we should have 
done differently, which should be part of a second 
generation simulation. First, we only recorded 
wave function information on a small number of 
time slices. We took a trial function iptrial(f) and 
recorded the C(i) resulting from Eqn. [|. Had we 
possessed wave function information on all sites, 
we could have reconstructed the wave function 
from the large t C(r, t) and used it as the tptrial ( r ) 
in Eqn. |4J. This might have given a better signal. 
Second, we used sources which were angular mo- 
mentum eigenstates only for the 3 P 2l 1 Pi, 3 £>3, 
and 1 Z?2 mesons, and the highest angular momen- 
tum baryons. The other states are merely eigen- 
states of m j . This was done to limit the number 
of propagators constructed. We could have fully 
reconstructed the sink angular momentum wave 



function, obtaining correlators which were eigen- 
functions of j and mj. We used nonrelativistic 
source and sink wave functions (upper Dirac com- 
ponents in Bjorken-Drell basis). This is probably 
good enough for heavier masses but might not 
be the optimal choice for very light quarks. Fi- 
nally, in excited state baryon spectroscopy one 
must deal with the fact that many states with 
the same quantum numbers are present in a mul- 
tiplct. For example, there are two j = 1/2 nucle- 
ons in the L — 1 [70] of SU(6). Separating these 
states presents a technical challenge. 

Our simulations used Wilson fermions on 16 4 
lattices at a coupling j3 — 6.0, in quenched ap- 
proximation, with 80 lattices for P waves and 50 
for D waves. They were performed at the Pitts- 
burgh Supercomputing Center. 

We present a picture of a D-wave hadron in Fig. 
[j]; a slice through a 3 D% meson. The light regions 
are where the wave function is large and positive; 
in the dark regions it is large and negative. The 
characteristic quadrupole lobes are obvious. The 
orbitally excited mesons and baryons are much 
larger than the S-wave states. 



Figure 1. Portrait of a 3 D^ meson. 



We were able to see a hint of P wave fine struc- 
ture splitting at our heaviest quark mass. It is 
shown in Fig. |[ With a nominal lattice spacing 
oil/ a — 2 GeV there is a remarkable similarity to 
charmonium, where the 3 Pq state is at 3415 MeV 



and the other states are at about 3500 MeV. 



Figure 2. Fine structure splitting at k = .1300. 



As a global way of presenting our results we 
show the "Wilson fermion wallet card" in Fig. pi 
We show S-wave mesons and baryons along with 
the 3 P 2 , N(5/2), 3 D 3 , and N(7/2) excited states 
(labelled "p'\ "P", "d", and "D".) 



lattice spacing by fitting the extrapolated masses 
to the V(3095) and x(3555) masses. This gave 
a charm hopping parameter of K = .1224 and 
an inverse lattice spacing of 1/a — 1790 MeV. 
(Determining the lattice spacing from our proton 
mass would give 1/a = 1710 MeV; from the rho, 
1/a = 2264 MeV.) The extrapolated common D- 
meson mass is then 3.99(16) MeV, where the error 
is only from the extrapolation. The 3 Di cc state 
is at 3.77 GeV but its mass is influenced by the 
nearby DD threshold. Model calculations H of 
D-wave states (some of which are narrow since 
their decays to DD are forbidden) give masses of 
3.81 to 3.84 GeV. At this value of the lattice spac- 
ing our 3 P2 — 3 Po mass splitting at k = .1300 is 63 
MeV; in charmonium the corresponding number 
is 145 MeV. 

This is a little different from the Fermilab pro- 
cedure of extracting a lattice spacing from the l P\ 
splitting from the center of gravity of the S-wave 
mesons. We made our choice because we have a 
better signal in the 3 P2 channel; anyway, the most 
robust experimental P-wave spectroscopy pretty 



Figure 3. Wilson spectroscopy vs. 



We extracted a lattice spacing by lattice deter- 
minations of the 3 Pi and 3 5i states at k — 0.130, 
0.145, and 0.152 and extrapolating their masses 
linearly in n. (We used the data of Ref. [|| for 
the n = .1300 vector meson). We determined the 



3. PHENOMENOLOGY 

While wave functions in a smooth gauge have 
clearly demonstrated their worth as interpolating 
fields for spectroscopy, it is still not clear if they 
have any real physical significance. We compared 
charge radii computed using wave functions with 
form factor data. The calculation is simple: for 
mesons one computes 



<r 2 > = E e > 



d 3 r#(r)| 2 (r/2) 5 



(0) 



and there is an analogous equation for baryons. 



We found that 



computed from Eqn. }B was 



about a factor of four to six smaller than the ob- 
served pion or nucleon charge radius. However, at 
the same time, we saw that the ratio of neutron 
to proton charge radii was negative and about the 
same size as the experimental ratio. Apparently 
this negative neutron charge radius arises because 
the two d quarks have a larger relative separation 
than the u quark does from either of the d's. This 
effect is easy to explain as being due to the color 



hyperfine interaction (compare Ref. 
show an example of this ratio in Fig. 



iF 



m). We 



Figure 4. Neutron/proton charge radius ratio. 



4. CONCLUSIONS 

The methods we describe here can and should 
be used for a large scale calculation of the spec- 
troscopy of orbital excitations. We believe our 
signal for P-wave states is only limited by statis- 
tics. However, the physical meaning of the wave 
function remains obscure. 
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